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We show that nonhnear response of a quantum oscillator displays antiresonant dips and resonant 
peaks with varying frequency of the driving field. The effect is a consequence of special symmetry 
and is related to resonant multiphoton mixing of several pairs of oscillator states at a time. We 
discuss the possibility to observe the antiresonance and the associated multiphoton Rabi oscillations 
in Josephson junctions. 
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Many resonant nonlinear phenomena are described by 
the model of a nonhnear osciUator in a resonant field. Ex- 
amples include collisionless dissociation of molecules IJ , 
dispersive optical bistability cyclotron resonance of 
a relativistic electron Q , resonant enhancement of acti- 
vated escape and tunnehng Q in Josephson junctions, 
and recently discovered hysteresis in Josephson junctions 
15 and nanomechanical resonators • 

A weakly nonlinear oscillator is a multi-level quantum 
system with nearly equidistant energy levels E^. There- 
fore a periodic force of frequency lof can be nearly res- 
onant for many transitions at a time, i.e., fiLOp can be 
close to the interlevel distance — £'„ for many n. 

This makes an oscillator convenient for studying multi- 
photon Rabi oscillations. They arise when the spacing 
between remote energy levels n and m coincides with 
the energy oi n — m photons. En — Em — (n — m)hLOp 
0. The multiphoton transition amplitude is resonantly 
enhanced, because the m n transition occurs via a 
sequence of virtual field-induced transitions k k + 1 
(with TO < fc < n — 1), all of which are almost resonant. 

In this paper we show that multiphoton transitions 
in the oscillator are accompanied by a new effect, an 
antiresonance of the response. When the frequency of 
the driving field adiabatically passes through a resonant 
value, the vibration amplitude displays a sharp minimum 
or maximum, depending on the initial conditions. We 
argue that the antiresonance and the multiphoton Rabi 
oscillations can be observed in such macroscopic systems 
as Josephson junctions and nanomechanical resonators. 

The multiphoton antiresonance is a consequence of two 
interesting properties of a driven oscillator. First, for res- 
onant photon frequencies, simultaneously, not one but 
several pairs of states turn out to be in resonance. Sec- 
ond, the amplitudes of forced vibrations in the resonating 
states coincide with each other, in the neglect of multi- 
photon mixing. For not too strong fields, these ampli- 
tudes are determined by nonresonant field-induced cou- 
pling of neighboring Fock states of the oscillator. The res- 
onant multiphoton mixing leads to level splitting, which 
strongly differs for different resonating pairs. It is this 
difference that results in the dips (peaks) in the vibration 
amplitude as adiabatically passes through resonances. 

In the semiclassical picture, resonant multiphoton 



transitions correspond to tunneling between Floquet 
states of the oscillator with equal quasienergies [the 
quasienergy e gives the change of the wave function ip{t) 
when time is incremented by the modulation period tf, 
^]{t+TF) = exp(— j£TF/?i)V'(i)]. The occurrence of equal- 
quasienergy states is related to the bistability of forced 
vibrations of a classical oscillator 8] . There is similarity 
between the oscillator tunneling and the tunneling of a 
particle in a static double-well potential, with the poten- 
tial minima being analogs of the stable states of forced 
vibrations. However, the latter states are not separated 
by a static barrier. Tunneling of a driven oscillator [p is 
a carefully studied example of dynamical tunneling |lOj. 
As we show, the WKB analysis gives an important in- 
sight into the origin of the antiresonance. 

The Hamiltonian of a driven nonlinear oscillator with 
mass M — 1 has the form 

H{t) - ^ 1^2g2 ^ 1^^4 _ qAcos{uJFt). (1) 

We assume that the driving field is nearly resonant, i.e., 
the frequency detuning 5lo is small, 

\5ijj\<^LjJf, 6uj = ujf ~ ojq- (2) 

We consider not too large amplitudes of the driving 
field A, so that the oscillator anharmonicity is small, and 
in particular \'y\q^ ^ loq for typical q. We also assume 
that 7 and 5uj have the same sign. This is required for the 
amplitude of classical forced vibrations to display hys- 
teresis as a function of the force amplitude A. If there 
is a cubic term ag'^/S in the potential energy of the os- 
cillator, its major effect of interest for this paper is the 
renormalization 7— i-7 — ^(a/wo)^ 0- 

To study quantum dynamics, we will write the Hamil- 
tonian in terms of the raising and lowering operators 
of the oscillator , a, and switch to the rotating frame 
with a canonical transformation U{t) — exp(— zwf a^at). 
The transformed Hamiltonian Hq = W{t)H{t)U{t) - 
ihU^ {t)U (t) is time-independent in the rotating wave ap- 
proximation (RWA), 

Hq = —5uon -f —Vfi{n + 1) — / (a + a^) , fi — a^a, 
V = 3h^/ALol f = {Stuuoy'^^A. (3) 
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The terms oc V, / that contain fast-oscillating factors oc 
exp{zLikojpt) with k — 2,4 were disregarded in Eq. 
In the expression for Hq and in what follows h — 1. 

The eigenvalues of the Hamiltonian Q £„ give the 
quasienergies of the driven oscillator. In the limit of weak 
driving their spectrum is particularly simple, 

Sn = -nSuj + Vn{n + 1) /2 (/ ^ 0). (4) 

We will study multiphoton resonance for the ground 
state of the oscillator, En — Eq = Nujp, or equivalently 
Eo = £n- From Eq. for small / and given N the 
resonance occurs for 

6uj = SujN = V^(iV + l)/2. 

Remarkably, at resonance all quasienergy levels Q with 
n < N are pairwise-degenerate, En-u = £n E3- Equiva- 
lently, Em-tl — En = {N — 2n)ujF- One can see that the 
degeneracy is not lifted by the lowest-order (c>c /^) field- 
induced level shift [except for the levels n — {N zt l)/2 
for odd N and n = (iV/2) ± 1 for even N]. As shown 
below, it persists for all / in the WKB approximation in 
the neglect of tunneling. 

The response of the oscillator to the field is character- 
ized by the expectation value of its coordinate q. If the 
oscillator is in an eigenstate \n) of the Hamiltonian J^Jl, 

g„ = (2wo)"^/^a„e-''^'=-* + c.c, a„ = {n\a\n). (5) 

To first order in the field, the reduced amplitude of 
forced vibrations a„ is 

a„ = -fSuj/{{Suj - Vn)[Suj - V{n + 1)]). (6) 

At multiphoton resonance, where Slu — 6ujn, the vibra- 
tion amplitudes in the resonating states coincide with 
each other, aAr_„ = a„ for < n < N/2. 

Multiphoton mixing leads to splitting of the 
quasienergy levels and the vibration amplitudes. It can 
be calculated by diagonalizing the Hamiltonian (PJ and 
is shown in Fig. ^ as a function of frequency detuning Soj. 
One of the involved resonating states is the ground state 
of the oscillator n = in the limit / — > 0. 

The minimal splitting of the levels Sq and sn is given 
by the multiphoton Rabi frequency Qn. For weak field 
it can be obtained from Eq. ||2J) by perturbation theory 
0. To the lowest order in f/SiVN 

nn^2f\2f/v\''-'N^my^^\ (?) 

For ^ 1 this expression becomes 

nn^V{f/M^N''/\2TTr^/\ 

/jv = |F|iV3/2exp(-3/2)/2. (8) 

The Rabi frequency depends on A^ exponentially, flji oc 
In the case A = 5 it is shown in Fig. 12 . One can 
see from this figure that Eq. ([7|) works well in the whole 
range of the field amplitudes ///at ^ 0.5. For larger 




FIG. 1: (a) Anticrossing of the quasienergy levels eo and es 
with varying frequency Suj for the field amplitude f/fs — 0.75, 
where the scaling field /s is given by Eq. (|HJ with N — 5. In 
the limit f the levels cross for Suj = Suis = 3V; (b) The 
reduced vibration amplitudes ao and as for the same adiabatic 
states. For (^i^s — Suj)/nR ^ 1 the states described by lines 
1 and 2 are close to the Fock states of the oscillator |0) and 
|5), respectively. 



fields r2/j depends on / much weaker than the asymptotic 
expression ® 0. 

The most interesting feature of Fig. ^ is the antireso- 
nant splitting of the amplitudes. It occurs at the adia- 
batic passage of 6uj through resonance, where the sys- 
tem switches between the ground and excited states. 
In particular, the amplitude displays an antiresonant 
dip if the oscillator is mostly in the ground state for 
{Su!—Su!n)/V< 1 or in the state A^ for {6llj—6lljn)/V > 1. 
The magnitude and sharpness of the dip are determined 
by ^Ir/V and depend very strongly on the field and A^. 
With decreasing ^Ir/V the dip (and peak) start looking 
like cusps located at resonant frequency. The amplitude 
splitting as function of the field is shown in Fig. |21 

To explain this behavior we note that, for Slo = Slon, 
the field leads to two major effects. One is resonant mix- 
ing of the wave functions into symmetric and antisym- 
metric combinations \n, N—n)± — (|n)o ± |A^ — rt)o) /V2 
with quasienergies e„± (|-)o are eigcnfunctions of the op- 
erator n — a^a). The second effect is nonresonant mixing 
of the states |n, N — n)± with different n, which leads to 
nonzero expectation values of the vibration amplitudes. 

To first order in /, the vibration amplitudes ao± = 
±(0, A^|a|0, A)± are determined by nonresonant mixing 
of the states |0, A^)± with \l,N - 1)± and \N + l)o. For 
comparatively weak fields, the level splitting ei+ — £i_ oc 
V,r{Suj/ f)"^ largely exceeds the splitting £o+ ^ ^o- = ^n- 
Then from perturbation theory ao+ — ag- oc /[(eq ^ 
- (eo - £i-)~^] oc {f/V)^-\ This scaling de- 
scribes the resonant small-field amplitude splitting in 
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FIG. 2: Upper panel: field-induced splitting of the 
quasienergy levels n = Q and n = N = 5 for resonant driving 
frequency, = Su!^ = 3V . The splitting gives the 5- photon 
Rabi frequency Qr. The dashed line shows the weak-field per- 
turbation theory ||7|l. Lower panel: splitting of the reduced 
amplitudes of forced vibrations in the corresponding Floquet 
states. The curve labelling coincides with that in Fig. 

Fig. 121 extremely well (the prefactor is determined by the 
admixture of states \n, N — n)± vi^ith n > 1 and will be 
discussed elsewhere). 

The simultaneous degeneracy of quasienergies and vi- 
bration amplitudes for many pairs of states in a broad 
field range can be shown analytically in the case where 
the oscillator dynamics is described by the WKB approx- 
imation. This approximation applies for 

A<1, X^V/{2Suj). (9) 

It is convenient to introduce the reduced coordinate 
and momentum of the oscillator in the rotating frame 

Q = {V/4SLoy/^{a + a^), P = -i{V/4SLo)^/^{a ~ a^) 

with the commutator [P, Q\ — ~iX. In these variables, in 
the neglect of terms (x A, the Hamiltonian ^ becomes 
Ho = 2{6u!)^V-'^[g{Q,P) - 1/4], where 

5(Q,P)-(Q' + P'-l)V4-/3'/'Q. (10) 

Here (3 = f'^V/{6uj)^ is the reduced field intensity. 

The function g{Q, P) is illustrated in Fig.|3 5uj g{Q, P) 
is the classical Hamiltonian in the RWA, it gives the 
quasienergy of the oscillator JJJ; Q,P are the canon- 
ical variables. The minimum and local maximum of 
g{Q,P) correspond to the stable states of forced vibra- 
tions. They coexist for < (3 < 4/27. For such (3, in 
a certain range of g there are two Hamiltonian trajecto- 
ries with the same g, one on the internal "dome" and the 



other on the external part of the surface g(Q,P). We call 
them, respectively, internal and external trajectories. 




FIG. 3: The classical quasienergy of the oscillator HUH . The 
plot refers to the reduced field /3 — 2/27. 

The external trajectory for given g{Q, P) = g has the 
form Q{t) = /3-i/2[X2(t) - g], with 

x{t) = <^^i-^-<^^)-;^^i<^^-^2Wu^ 

ci - C3 - (ci - C2)sn^u 

Here, sn u is the Jacobi elliptic function; the elliptic mod- 
ulus is m = (ci - C2)(c3 - Ci)/{ci - c^){c2 - C4), and 
u = [(ci— C3)(c2— C4)]^/^(5(jj t/2 is the appropriately scaled 
time. The coefficients ci > C2 > C3 > C4 are the roots 
of the polynomial /3(1 -I- 2x) — {x^ — gY {X = ci, C2 and 
X = C3, C4 gives the turning points Q = on the external 
and internal trajectories, respectively). 

The internal trajectory Q{t) is given by Eq. Hill) with 
u u + K + iK', where K = K{m) is the complete 
elliptic integral, and K' = K{1 — m). 

An immediate consequence of the analytical interrela- 
tion between the external and internal trajectories is that 
the periods of motion along them coincide The vi- 
bration frequency isuj{g) — TT[{ci—C'i){c2 — C4)Y^''^5u}/2K. 
When motion is quantized, Lo{g) gives the distance be- 
tween the energy levels. Therefore if, for some Su and (3, 
two levels that correspond to the external and internal 
trajectories coincide with each other, many levels will co- 
incide pairwise as well. Resonant multiphoton splitting 
of the levels is due to tunneling between the external and 
internal parts of the surface g{Q,P). 

In the WKB approximation the expectation value a„ 
in a quantum state |n) is given by the period-averaged 
value {Q)g of the coordinate Q on the appropriate classi- 
cal trajectory. The values of {Q)g turn out to be the same 
on the internal and external trajectories with the same 
g. Indeed, their difference can be written as a contour 
integral of in the u-plane. The contour is a parallelo- 
gram with vertices at u = 0, 2K, 3K + iK', and K + iK' . 
The function X'^ has one second order pole in this par- 
allelogram. The residue is cx ci -I- C2 + C3 -|- C4 = 0. This 
proves that the vibration amplitudes of the states with 
equal quasienergies are equal, in the neglect of tunneling. 

In order to observe the coherent multiphoton quantum 
effects, the Rabi frequency flu should exceed the deco- 
herence rate. In the RWA, relaxation of an oscillator can 
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often be described by the master equation for the density 
matrix p. In many cases of interest the major source of 
decoherence is noise or a quantum field that modulates 
the oscillator energy. Then the coupling operator is oc n. 
If the modulation is fast compared to 6uj, V , the master 
equation takes the form 

P^i[p,Ho]~t^p, fp^T^[h,[n,p]], (12) 

where characterizes the noise intensity. It gives the 
phase diffusion coefficient of the oscillator. 

From Eq. 112|l . for resonant Slu = (5a; jy the major effect 
of phase diffusion is decay of the Rabi oscillations. For 
ilji/SujN <C 1 the decay rate is F^A^^: this is simply the 
diffusion coefficient of the phase difference of the Fock 
states of the oscillator |0) and \N). More formally, this is 

the decay rate of the difference p++ — p , where p^i,' — 

,{0,N\p\0,N),, (z/,:/' = ±). 

The decoherence rate due to phase diffusion quickly in- 
creases with N. Our results demonstrate that the strong 
resonant amplitude splitting occurs already for N = 5. 
We emphasize that this is a coherent quantum effect. It 
is qualitatively different from the nonmonotonic field de- 
pendence of the stationary amplitude of a driven damped 
oscillator for nonzero temperatures [T3. [Tsl . 

The antiresonances in the vibration amplitude can be 
directly observed in Josephson junctions (JJ's), which are 
well described by the model of a nonlinear oscillator . 
The amplitude of forced vibrations of a JJ was recently 
measured and the vibration bistability was demonstrated 
1^ . The eigenfrequency and nonlinearity of a J J are con- 
trolled by the dc bias current /dc ■ When I^c is close to the 
critical current /q, i.e., rj ^ {Iq — Idc)/Io ^ li hi Eq. Q 
Wo = ujp{2r]y^'^ and V = —bhujp/iSEjrj, where ujp is 
the J J plasma frequency and Ej = hIo/2e. The scaling 
RF current In for a resonant A/'-photon transition, which 
corresponds to the scaling field /at in Eq. is 



In = (5/48)/cexp(-3/2)(2iVncc;p/i;j)3/2(277)-^/8. 

For the RF current /rf ^ In at resonant frequency 
ujp = ojQ + V{N + l)/2, the splitting of the vibration 
amphtude is strong, see Figs. Q El Roughly, In is in- 
versely proportional to Iq — /dc, as is also the distance 
V/2 between the resonant values LUp. 

Multiphoton Rabi oscillations in multilevel JJs can be 
also studied by measuring the rate of tunneling decay in- 
duced by resonant RF pulses of different length ^rf- The 
effective potential of a JJ is metastable, uJ^q^ (2 — aq^ /Z. 
The probability of tunneling escape from an excited state 
is much larger than from the ground state, and the higher 
the energy level the stronger is the difference. Therefore 
at multiphoton resonance the escape rate should oscillate 
as sin^(rifl,iR,F/2). This approach has been used to detect 
single-photon Rabi oscillations in strongly nonlinear JJ's 
that have a small number of metastable states 0| . 

We note that the quantum effects that we discuss can 
be studied also in nanomechanical systems where classi- 
cal bistability of forced vibrations was observed JJ- 

In this paper we have shown that multiphoton response 
of a quantum oscillator may display antiresonant dips 
(peaks) as a function of frequency. The effect is due to 
level anticrossing and coherent mixing of many oscillator 
states at a time. The shape and magnitude of the dips 
(peaks) strongly depend on the field. We find the mul- 
tiphoton Rabi frequency and discuss the possibility to 
observe the antiresonant response and multiphoton Rabi 
oscillations in Josephson junctions and nanomechanical 
resonators. 
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